Insights into functions of the H channel of cytochrome c oxidase from atomistic molecular dynamics simulations.

Proton pumping A-type cytochrome c oxidase (CcO) terminates the respiratory chains of mitochondria and many bacteria. Three possible proton transfer pathways (D, K, and H channels) have been identified based on structural, functional, and mutational data. Whereas the D channel provides the route for all pumped protons in bacterial A-type CcOs, studies of bovine mitochondrial CcO have led to suggestions that its H channel instead provides this route. Here, we have studied H-channel function by performing atomistic molecular dynamics simulations on the entire, as well as core, structure of bovine CcO in a lipid-solvent environment. The majority of residues in the H channel do not undergo large conformational fluctuations. Its upper and middle regions have adequate hydration and H-bonding residues to form potential proton-conducting channels, and Asp51 exhibits conformational fluctuations that have been observed crystallographically. In contrast, throughout the simulations, we do not observe transient water networks that could support proton transfer from the N phase toward heme a via neutral His413, regardless of a labile H bond between Ser382 and the hydroxyethylfarnesyl group of heme a In fact, the region around His413 only became sufficiently hydrated when His413 was fixed in its protonated imidazolium state, but its calculated pKa is too low for this to provide the means to create a proton transfer pathway. Our simulations show that the electric dipole moment of residues around heme a changes with the redox state, hence suggesting that the H channel could play a more general role as a dielectric well.

Core level binding energies of functionalized and defective graphene.

X-ray photoelectron spectroscopy (XPS) is a widely used tool for studying the chemical composition of materials and it is a standard technique in surface science and technology. XPS is particularly useful for characterizing nanostructures such as carbon nanomaterials due to their reduced dimensionality. In order to assign the measured binding energies to specific bonding environments, reference energy values need to be known. Experimental measurements of the core level signals of the elements present in novel materials such as graphene have often been compared to values measured for molecules, or calculated for finite clusters. Here we have calculated core level binding energies for variously functionalized or defected graphene by delta Kohn-Sham total energy differences in the real-space grid-based projector-augmented wave density functional theory code (GPAW). To accurately model extended systems, we applied periodic boundary conditions in large unit cells to avoid computational artifacts. In select cases, we compared the results to all-electron calculations using an ab initio molecular simulations (FHI-aims) code. We calculated the carbon and oxygen 1s core level binding energies for oxygen and hydrogen functionalities such as graphane-like hydrogenation, and epoxide, hydroxide and carboxylic functional groups. In all cases, we considered binding energy contributions arising from carbon atoms up to the third nearest neighbor from the functional group, and plotted C 1s line shapes by using experimentally realistic broadenings. Furthermore, we simulated the simplest atomic defects, namely single and double vacancies and the Stone-Thrower-Wales defect. Finally, we studied modifications of a reactive single vacancy with O and H functionalities, and compared the calculated values to data found in the literature.

Reduction potentials and acidity constants of Mn superoxide dismutase calculated by QM/MM free-energy methods.

We used two theoretical methods to estimate reduction potentials and acidity constants in Mn superoxide dismutase (MnSOD), namely combined quantum mechanical and molecular mechanics (QM/MM) thermodynamic cycle perturbation (QTCP) and the QM/MM-PBSA approach. In the latter, QM/MM energies are combined with continuum solvation energies calculated by solving the Poisson-Boltzmann equation (PB) or by the generalised Born approach (GB) and non-polar solvation energies calculated from the solvent-exposed surface area. We show that using the QTCP method, we can obtain accurate and precise estimates of the proton-coupled reduction potential for MnSOD, 0.30±0.01 V, which compares favourably with experimental estimates of 0.26-0.40 V. However, the calculated potentials depend strongly on the DFT functional used: The B3LYP functional gives 0.6 V more positive potentials than the PBE functional. The QM/MM-PBSA approach leads to somewhat too high reduction potentials for the coupled reaction and the results depend on the solvation model used. For reactions involving a change in the net charge of the metal site, the corresponding results differ by up to 1.3 V or 24 pK(a) units, rendering the QM/MM-PBSA method useless to determine absolute potentials. However, it may still be useful to estimate relative shifts, although the QTCP method is expected to be more accurate.

Influence of the [2Fe]H subcluster environment on the properties of key intermediates in the catalytic cycle of [FeFe] hydrogenases: hints for the rational design of synthetic catalysts.

Nature's recipe: A theoretical study analyzes how the environment of the [FeFe] hydrogenase's catalytic cofactor affects its chemical properties, particularly the relative stability of complexes with bridging and terminal hydride ligands (see picture; Fe teal, S yellow, C green, N blue, O red, H gray). The results help to elucidate key rules for the design of bioinspired synthetic catalysts for H(2) production.

QM/MM-PBSA method to estimate free energies for reactions in proteins.

We have developed a method to estimate free energies of reactions in proteins, called QM/MM-PBSA. It estimates the internal energy of the reactive site by quantum mechanical (QM) calculations, whereas bonded, electrostatic, and van der Waals interactions with the surrounding protein are calculated at the molecular mechanics (MM) level. The electrostatic part of the solvation energy of the reactant and the product is estimated by solving the Poisson-Boltzmann (PB) equation, and the nonpolar part of the solvation energy is estimated from the change in solvent-accessible surface area (SA). Finally, the change in entropy is estimated from the vibrational frequencies. We test this method for five proton-transfer reactions in the active sites of [Ni,Fe] hydrogenase and copper nitrite reductase. We show that QM/MM-PBSA reproduces the results of a strict QM/MM free-energy perturbation method with a mean absolute deviation (MAD) of 8-10 kJ/mol if snapshots from molecular dynamics simulations are used and 4-14 kJ/mol if a single QM/MM structure is used. This is appreciably better than the original QM/MM results or if the QM energies are supplemented with a point-charge model, a self-consistent reaction field, or a PB model of the protein and the solvent, which give MADs of 22-36 kJ/mol for the same test set.

Proton Transfer at Metal Sites in Proteins Studied by Quantum Mechanical Free-Energy Perturbations.

Catalytic metal sites in enzymes frequently have second-sphere carboxylate groups that neutralize the charge of the site and share protons with first-sphere ligands. This gives rise to an ambiguity concerning the position of this proton, which has turned out to be hard to settle with experimental, as well as theoretical, methods. We study three such proton-transfer reactions in two proteins and show that, in [Ni,Fe] hydrogenase, the bridging Cys-546 ligand is deprotonated by His-79, whereas in oxidized copper nitrite reductase, the His-100 ligand is neutral and the copper-bound water molecule is deprotonated by Asp-98. We show that these reactions strongly depend on the electrostatic interactions with the surrounding protein and solvent, because there is a large change in the dipole moment of the active site (2-6 D). Neither vacuum quantum mechanical (QM) calculations with large models, a continuum solvent, or a Poisson-Boltzmann treatment of the surroundings, nor combined QM and molecular mechanics (QM/MM) optimizations give reliable estimates of the proton-transfer energies (mean absolute deviations of over 20 kJ/mol). Instead, QM/MM free-energy perturbations are needed to obtain reliable estimates of the reaction energies. These calculations also indicate what interactions and residues are important for the energy, showing how the quantum system may be systematically enlarged. With such a procedure, results with an uncertainty of ∼10 kJ/mol can be obtained, provided that a proper QM method is used.

Calculated reaction cycle of cytochrome c oxidase.

The catalytic cycle of cytochrome c oxidase has been simulated by means of quantum mechanical calculations. The experimental energetics of the catalytic cycle is nearly reproduced. The atomic structures of the intermediates are suggested. In particular, the structures of nonactive "resting" intermediates are proposed.